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ABSTRACT 


Theoretical  Rayleigh  wave  dispersion  curves 
for  three  different  types  of  sub-oceanic  basement 
layering  are  presented.  Previous  conclusions 
concerning  the  dispersion  of  Rayleigh  waves  across 
ocean  basins  are  re-examined  in  the  light  of  the 


new  data. 
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I.  Introduction 

Recent  Investigations^-  have  shown  that  the 
characteristics  of  propagation  of  Rayleigh  waves 
across  ocean  basins  can  be  explained  by  the  effect 
of  a  layer  of  water  overlying  ultrabasic  rock  having 
an  equivalent  compressional  velocity  of  slightly 
under  7,90  km/sec.  In  this  earlier  work  the 
theoretical  problem  was  reduced  to  its  simplest  form 
by  freacing  the  water  and  sediment  as  a  single 
liquid  layer  and  the  basement  rock  as  a  single 
homogeneous  layer  of  infinite  thickness. 

Since  the  theoretical  dispersion  curve  from  this 
idealized  structure  agreed  remarkably  well  with  the 
observed  dispersion  curve  it  was  concluded  that  any 
layering  present  in  the  basement  rocks  involved  only 
very  moderate  contrasts  In  elastic  properties.  It 
was  pointed  out  that  the  small  dispersion  observed 
for  Love  waves  along  ocean  paths  demanded  a  super¬ 
ficial  layer  in  the  basement  rocks,  but  that  the 
elastic  properties  of  the  layer  should  be  much  closer 
to  those  of  the  substratum  than  is  the  case  for  the 
continents.  Seismic  refraction  measurements 
the  Atlantic  and  Pacific  ocean  basins  have  given  the 
thiclmess  and  elastic  constants  of  the  basement  layer 
in  a  number  of  places.  The  layer  is  roughly  as  thick 
as  the  ocean  is  deep  and  has  a  velocity  for 
compressional  waves  varying  between  6.  1/4-7  km/sec. 
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usually  overlying  more  basic  rock  with  a  velocity 
range  of  about  7. 9-8. 2  km/sec. 

In  the  present  paper  calculations  of  the  effect 
of  the  basement  layering  on  Rayleigh  wave  dispersion 
curves  are  presented  for  three  different  types  of 
basement  layering: 

Case  I  <%o=  1.52  km/sec,  =  5.5  km/sec,  =  8.1 

km/sec,  H>  =  Hc,  H %  -  oo  ,  -  2.67  pQ  ,  pz  =  3.0  ; 

Case  II  Oio-  1.52  km/sec.  oCy  ■  7.90  km/sec,  H/  -CO  . 
p/  =3.0  p0\  Case  III  OiQ  =  1.52  km/sec  ; 

0C[  =  6.9  km/sec,  =  8.1  km/sec,  K,  =  He,  Fg  =  CO  , 
=  2.67  =  3.0  .  The  subscripts  0,1,2 

refer  respectively  to  the  liquid  layer,  first  and 
second  basement  layers,  and  Poisson1 2 * 4 s  constant  is 
taken  as  1/4. 

Comparison  of  the  new  theoretical  curves  with 
observed  dispersion  is  made.  Applications  to 
microseism  propagation  will  be  discussed  in  a  future 
paper. 


1.  Maurice  Swing  and  Frank  Press  "Crustal  Structure  and 
Surface  Wave  Dispersion  Part  II,  Solomon  Islands  Earth¬ 
quake  of  29  July  1950" ,  in  press  Bull.  Seism.  Soc.  Amer. 

2. C.  B.  Officer,  Jr.,  Maurice  Ewing  and  Paul  Wuenschel, 
"Seismic  Refraction  Measurements  in  the  Atlantic  Ocean 
Basin.  Part  IV,  Bermuda,  Bermuda  Rise  and  Wares  Deep, 
in  press  Bull.  Geol.  Soc.  Amer. 

3- Maurice  Ewing,  George  Sutton,  C.B.  Officer,  Jr., 
"Seismic  Refraction  Measurements  in  the  Atlantic  Ocean 
Basin,  Part  VI-  Typical  Deep  Stations  in  the  North 
Atlantic  Basin,  in  preparation. 

4.  Russel  W.  Raitt,  "Seismic  Refraction  Studies  of  the 
Pacific  Ocean  Basin",  paper  presented  before  annual 
meeting  Seism.  Soc.  America,  March  24,  1951. 
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II  Derivation  of  Period  Equation 


In  order  to  find  the  phase  and  group  velocity 

in  terms  of  the  parameter  kH  (wave  number  x  the  depth 

of  the  liquid  layer)  for  a  system  consisting  of  a 

liquid  superposed  on  two  solid  layers,  we  consider  the 

propagation  of  a  train  of  plane  waves  of  the  RayLeigh 

type.  We  make  use  of  notations  represented  in  Figure 

1,  where  are  the  densities,  and  ti  elastic 

constants,  oc .  ~ and  the  veloci- 

1  r  fc  ¥  y< 

ties  of  compress ional  and  shear  waves  respectively. 

As  usual  the  horizontal  and  vertical  components 
of  the  displacement,  u  and  w,  are  represented  by  the 
formulae 


dz 


and  the  functions  and  ^  should  satisfy  the  wave 
equation,  i.e. 


(W^j r-o,  (v >+%)?*(/  m 

The  stresses  are  given  by  the  equations 

A, -I n  ,  fa -fife  -  M )  w 


Bearing  in  mind  the  assumed  layering  for  cases  I  and 
III,  w©  take  =  H;  =  H;  the  equation  of  the  free 
surface  is  z  =  -  H.  It  is  obvious  that  the  final 
results  of  this  investigation  may  be  easily  changed  In 


order  to  correspond  to  the  more  general  case  H, 

Now  the  boundary  conditions  of  the  problem  may 
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be  written  in  the  form 


Ps<> 

at 

z—H 

(Pzz), 

=  (p~ ) 

l  >  (P7*X~- 

■°  u 

z=0 

(0 

K  •*, 

(Pzz), 

,  (fzz  \  3 

PJzx)i  > 

z=H . 

(6) 

(  at 

We  have  8  boundary 

conditions 

and  the 

expressions 

for 

the  functions  and  j/;  will  also  contain  8  unknown 
coefficients.  Omitting  the  time  factor  we  write  these 
functions  in  the  form 


ft  ={d  ey  (-/.  z)  +  Cejy>  fa  zjJ  expf  Ux)  <  7> 

ft-jD&xpCft?)  +£  exf/l,  zjj  t/p  f-cWx)  (*) 
ft  *  [F ejcp(-fc  *)  exp  (-inx.)  M 

ft  =  {Mcxfijf,  z)+t(-  exp(-fo  2))  (/o) 

ft  ~  2)  «J <jf>(-c*K)  O') 


Upon  eliminating  the  constants  B,C,D,E,F,M,N,P  between 
the  eight  boundary  equations  (4) -(6)  we  can  obtain  the 
period  equation  relating  phase  velocity  to  wave  number. 
The  method  of  elimination  is  similar  to  that  described 
in  a  previous  paper  ®  in  which  the  liquid  layer  was 
taken  as  infinitely  thick. 


5  W.  S.  Jardetzky  and  Frank  Press,  ”Rayleigh  Wave 
Coupling  to  Atmospheric  Compressional  Waves”,  in  press. 
Bull.  Seism.  Soc.  Amer. 
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For  briefness  we  make  the  following  definitions 

C  i  /  C 

and 


V„s  k  > 


0%) 


n-.-'/M?  ,  '«,=^R5|  'W| 

v'(Rrsk  ’ 

Making  use  of  Lee^  6  notation  we  obtain 

y -v;$  <■%%-!)-'? 


(n) 


c(‘ 


7  =  A  -  /*s£  -  2(p  -/)=  X-V. 

^  ft  r"  /31  r°  A  A/  7  A 


Putting  (see  equations  34-36  of  the  previous  paper^) 

6i*-xz-*k»y  Vy7  6-3=^«v^-aY^ 

urU 

i=Jt(z-vS$)6-, 

<-<*-'&?$■  % \f; + ^ v-i  (,6) 


&  «vy'v£^  + 

/«  '  ^  3.  / 


fy£*m.ny 

*  e,  1  *  *  ^ 


4s«-f£)G,  -*<5* 


6  A.W.  Lee,  "The  effect  of  Geologic  Structure  upon 
Microseismic  Disturbance",  M*N«R.A*S*  (Geophys.  Suppl.) 
Vol.  3,  pp.  83-105,  1932 
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/'_  _  ^  Vs  -%  tanfL(n^h) 

<  "  p‘  n0n^  '•  },  K 


A 

4"- 

4"=-^r  ^fPfi'  tanktn.Kh) 


(17) 


the  period  equation  takes  the  form 


£  +  (f ' %onJ(n,Kf1)i*Hl  (njK.h)-(-  l^$ot)k(nl  xh)c&sh  (nh  rtk ) 

4  l' orth'nli£)!>!n.k(nittk)+£t/c.oih(nl/<!fi)c<K4.(n}K:H)  ,gs 

"^(n,  kff)s/nk(n}  *s<aA  (n,  iik)cM  (*}  *  «) 

+  t^otik  (n/KHjSifik(niick)+  Cj  /ckfj&ih  P4  t*nh  (njk) = 0 . 
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The  period  equation  for  the  more  general  case,  where 
the  depth  of  the  liquid  layer  differs  from  that  of 
the  upper  solid  layer,  will  have  the  same  form  as 
equation  ( /$),  the  factor  H  in  tanh  (n0kH)  being  re¬ 
placed  by  a  new  value,  say  E0 . 

It  may  be  noted  that  the  period  equation  for  the 
present  problem  can  be  obtained  from  that  for  an 
infinite  liquid  layer5  by  replacing  with  tanh(/?0  tfhj. 

The  numerical  procedure  for  computing  kH  as  a 
function  of  Vc  =  ^/<xo  consists  of  assuming  values  of 
and  obtaining  kH  from  {/S )  by  successive 
approximations.  Group  velocity  is  obtained  by  a 
graphical  differentiation  of  the  phase  velocity  curves. 
Periods  are  related  to  kH  through  the  equation  T  * 

~%irH /«(,*#  o/<xo  ■ 

The  derivation  of  the  period  equation  and  a 
description  of  methods  of  computation  for  case  II 

rj 

were  given  in  previous  papers.' 

Haskell7 8  recently  derived  in  matrix  form 
appropriate  equations  for  computing  Rayleigh  dispersion 
in  three  layers. 

7.  Frank  Press,  Maurice  Ewing,  and  Ivan  Tolstoy,  ”The 
Airy  Phase  of  Shallow  Focus  Submarine  Earthquakes”,  Bull. 
Seism.  Soc.  Vol.  40  pp  111-148,  1950. 

Q.  Norman  A.  Haskell,  flThe  Dispersion  of  Surface  Waves  on 
Multi-Layered  Media”,  Geophysical  Research  Paper  No.  9,  Air 
Force  Cambridge  Research  Center,  Cambridge,  Mass.,  1951. 


9 


III  Discussion 

The  numerical  results  for  phase  velocity  are 
presented  for  the  lowest  modes  of  cases  I,  II  and  III 
in  Tables  1  and  2.  Group  velocity  curves  for  these 
three  cases  are  presented  in  dimensionless  form  in 
Figure  2*  Group  velocity  curves  for  periods  between 
15  and  40  sec  are  presented  for  cases  I  and  III  in 
Figure  3  for  an  ocean  depth  H  =  5*57  km.  As  in 
previous  papers  we  have  taken  this  value  of  H  to  in¬ 
clude  the  thickness  of  unconsolidated  deep  sea  sediments. 
From  examination  of  these  curves  the  following 
conclusions  can  be  reached, 

1-  Gases  II  and  III  are  experimentally  in¬ 
distinguishable  from  each  other  in  the  range  15-40 
sec.  For  T  >  40  sec  the  two  curves  will  diverge. 

Previous  results ^  based  on  Rayleigh  wave  dispersion 
indicated  that  the  Pacific  Ocean  basin  is  underlain 
by  rOCK  having  an  equivalent  velocity  of  slightly 
under  7,90  km/sec.  It  is  now  seen  that  a  layered 

I 

basement  approximating  the  properties  assumed  for  case 
III  with  H  =  5,57' km  explains  the  observed  dispersion 
equally  well.  This  is  consistent  with  refraction 
measurements  in  the  Atlantic  and  Pacific  Oceans. 

2-  The  sensitivity  of  the  technique  of 
comparing  observed  Rayleigh  wave  dispersion  curves  with 
theoretical  curves  to  ascertain  suboceanic  basement 
structure  can  be  estimated  from  Figure  3,  Since 
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changes  in  H  simply  shift  the  group  velocity  curve 
to  the  right  or  left,  it  can  be  seen  that  these  curves 
are  more  sensitive  to  changes  in  velocity  of  the  first 
basement  layer  in  the  range  23  <  T  <40  sec*  For  a 
well  located  earthquake  having  a  predominatly  oceanic 
path2,  the  scatter  of  observed  Rayleigh  wave  velocities 
in  this  period  range  is  less  than  0*1  km/sec.  One  can 
conclude  that  differences  in  structure  of  the  first 
basement  layer  from  the  conditions  assumed  in  case  III 
with  H  =  5.57  km,  can  be  detected  by  analysis  of 
Rayleigh  wave  dispersion  providing  the  layer  thickness 
is  at  least  as  great  as  the  ocean  is  deep  when  the 
velocity  contrast  is  as  great  as  that  between  cases  I 
and  III* 

It  Is  to  be  hoped  that  further  investigations  will 
reveal  techniques  for  correcting  for  variations  in 
water  depth  along  paths,  refraction  at  the  continental 
border,  so  that  determination  of  crustal  structure  can 
be  made  over  most  oceanic  areas  from  the  records  of 
existing  seismological  stations* 

We  are  grateful  to  Professor  Maurice  Ewing  for 
many  helpful  discussions  on  this  problem* 


Table  I  -  Phase  Velocity  Values 


Case  I 

Case 

III 

c/x. 

kH 

£/<*<, 

kH 

1.00 

10.667 

1.00 

47.206 

1.05 

4.204 

1.05 

4.431 

1.075 

3.483 

1.075 

3.646 

1.10 

2.990 

1.10 

3.168 

1.20 

2.130 

1.20 

2.220 

1.30 

1.726 

1.30 

1.782 

1.40 

1.471 

1.40 

1.513 

1.50 

1.289 

1.50 

1.326 

1.60 

1.156 

1.60 

1.185 

1.65 

1.099 

1.65 

1.126 

1.70 

1.049 

1.70 

1.073 

1.75 

1.003 

1,75 

1.026 

1.80 

.961 

1.80 

.982 

1.85 

.921 

1.85 

.942 

1.90 

.883 

1.90 

.906 

1.95 

.850 

1.95 

.871 

2.00 

.818 

2.00 

.839 

.205 

.789 

2.05 

.808 

.210 

.759 

2.10 

.781 

.220 

.704 

2.20 

.728 

.230 

.656 

2.30 

.678 

.240 

.602 

2.40 

,629 

.250 

.537 

2.50 

.579 

.260 

.442 

2.60 

.510 

.270 

.283 

2.70 

.380 

.273 

.222 

2.73 

.323 

.275 

.177 

2.75 

.268 

.277 

.123 

2.77 

.204 

.280 

.061 

2.80 

.098 

Table  II.  Phase  Velocity  Values  Case  II 


c  Mo 

kH 

1.00 

32.863 

1.04 

5.0945 

1.08 

3.6337 

1.12 

2.96086 

1.16 

2.54963 

1.20 

2.26415 

1.24 

2.05061 

1.28 

1.88272 

1.32 

1.74606 

1.36 

1.63187 

1.40 

1.53451 

1.44 

1.45016 

1.48 

1.37609 

1.52 

1.31033 

1.56 

1.25141 

1.60 

1.19847 

1.64 

1 . 14975 

1.68 

1.10542 

1.72 

1.06461 

1.76 

1.02685 

1.80 

.99179 

1.84 

.95896 

1.88 

.92825 

c/<X0 

kH 

1.92 

.89932 

1.96 

.87201 

2.00 

.84611 

2.04 

. 82 144 

2.08 

.79786 

2.12 

.77521 

2.16 

.75337 

2.20 

.73218 

2.24 

.71149 

2.28 

.69114 

2.32 

.67094 

2.36 

.65065 

2.40 

.62998 

2.44 

.60850 

2.48 

.58559 

2.52 

.56029 

2.56 

.53097 

2.60 

.49463 

2.64 

.44532 

2.68 

.36854 

2.72 

.23092 

2.74 

.  12 196 
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Figure  i.  Nomenclature  and  coordinate  system. 


Figure  2.  Dimensionless  group  velocity  curve  for  cases  I,  II  and  III. 
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